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Abstract 

The majority of the highest energy cosmic rays are thought to be electrically 
charged: protons or nuclei. Charged particles experience angular deflections as they 
pass through galactic and extra-galactic magnetic fields. As a consequence correla- 
tion of cosmic ray arrival directions with potential sources has proved to be difficult. 
This situation is not helped by current data samples where the number of cosmic 
rays/source are typically < 0(1). Progress will be made when there are signifi- 
cantly larger data samples and perhaps with better catalogs of candidate sources. 
This paper reports a search for correlations between the RXTE catalog of nearby 
active galactic nuclei, AGNs, and the published list of ultra-high energy cosmic 
rays from the AG AS A experiment. Although no statistically significant correlations 
were found, two correlations were observed between AGASA events and the most 
inclusive category of RXTE AGNs. 
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1 Introduction 



Perhaps the primary goal of all experiments studying the highest energy cos- 
mic rays is to find the source of these particles. While circumstantial evidence 
may favor one type of source over another, demonstration of a clear correlation 
between the direction of cosmic rays and their sources is arguably essential. 
Unfortunately for electrically charged cosmic rays, galactic magnetic fields, 
and for the highest energy cosmic rays extra-galactic magnetic fields, cause 
angular deflections that can blur the correlation between cosmic ray arrival 
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direction and source direction. If the sources as viewed from the earth are ex- 
tended [1,2] the problem is even more difiicuh. Unless otherwise noted, for this 
paper we assume compact (point-like) sources for the highest energy cosmic 
rays. 

If the angular blurring from magnetic fields is small [3] {i.e. not significantly 
greater than the experimental angular resolution) and/or for neutral primaries, 
then experiments should observe cosmic rays that cluster in arrival direc- 
tion [4,5], and/or that correlate with potential astronomical {e.g. BL Lac) 
sources [6,7,8,9,10]. For nearby sources, where experiments should detect mul- 
tiple cosmic rays/source, event clusters provide bounds on the cosmic ray 
source density [11,12,13] potentially favoring one type of source for the high- 
est energy cosmic rays over another. However at this time the situation is 
less than clear as some results [14,15] question the significance of the reported 
clusters and/or some of the BL Lac correlations [10,16]. 

If defiections of charged cosmic rays by extra-galactic magnetic fields are not 
small [17], then lower energy, E, cosmic rays should experience the greatest 
angular deflections. Unfortunately small experiment data samples and a cos- 
mic ray flux oc have often caused studies to retain cosmic rays to ener- 
gies, Ethresh, well below GZK [18] energies [4]. Furthermore deflections of the 
highest energy cosmic rays even by our galactic magnetic field can be sub- 
stantial [19,20,21]. As magnetic defiections scale proportional to the charge of 
the primary cosmic ray, nuclei in the cosmic rays may have significant de- 
flections. Although most searches have looked for clustering and/or source 
correlations on small angular scales, studies at larger angular scales have also 
found evidence for clustering and/or source correlations [22,23,24]. Certainly 
the angular scale of cosmic ray clusters and the magnitude, and thus relevance, 
of the defiections of ultra-high energy cosmic rays by magnetic fields is not 
universally agreed to at this time. 

In the future, significantly larger data samples will allow analyses to increase 
Ethresh while retaining the number of observed cosmic rays/source (for nearby 
sources) > 0(1). However another possibility is to exploit catalogs of candi- 
date sources. With a catalog of source directions, cosmic rays can be effectively 
correlated with sources even when if magnetic field deflections are "not small" 
and/or when the the number of observed cosmic rays per source is < 1 allow- 
ing searches with existing data samples. That said, catalog based studies are 
limited by the completeness of the source catalog and the relevance (or not) of 
that class of astronomical source to the production of the highest energy cos- 
mic rays. Often conjectured astrophysical sources include gamma ray bursts, 
GRBs, and/or active galactic nuclei, AGNs[25]. 

This paper reports a search for correlations between a catalog of nearby 
AGNs [26] and the pubhshed list of ultra-high energy cosmic rays from AGASA [4] . 
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The components of our analysis are listed in Section 2. Issues that relate to 
data and AGN selection are given in Section 3. The cosmic ray-AGN compar- 
ison results are given in Section 4. Section 5 summarizes this study. 



2 Analysis Components 

Our comparison of ultra-high energy cosmic rays and a catalog of AGNs in- 
cludes three components: the RXTE catalog of AGNs, the AGASA list of 
cosmic rays, and a Monte Carlo sample of uniformly distributed cosmic rays 
generated to match the experimental acceptance of AGASA. 

The catalog of nearby AGNs [26] results from the Rossi X-ray Timing Explorer, 
RXTE, all-sky slew survey [27] sensitive to sources of hard X-rays (3-20keV). 
The survey excluded the galactic plane {\b\ > 10°) but covered ~ 90% of the 
remaining sky. X-ray sources were located to better than 1° and then correlated 
with known astronomical objects. The efficiency for AGN identification was 
estimated to be ~ 70% with somewhat higher efficiency for northern AGNs 
(~ 87%) and somewhat lower efficiency for southern AGNs (~ 60%) [26]. 
The resulting catalog provides source directions and probable source distances 
and intrinsic X-ray luminosities, I/3_2o- The catalog is best for nearby AGNs 
as RXTE signal thresholds significantly reduced the efficiency for detecting 
distant sources; additional details are given below. 

The list of ultra-high energy cosmic rays comes from published AGASA data 
[4]. 

The Monte Carlo sample of uniformly distributed cosmic rays was generated 
according to a cos(^)sin(^) distribution in local zenith angle, 9 < 45°, and 
uniform in local azimuth. Events were then transformed to celestial right as- 
cension and declination assuming constant detector aperture with time. 

Correlations, between the AGASA events and the catalog of AGNs from 
RXTE, would appear as an excess at small angular separations in compar- 
ison to the Monte Carlo sample of simulated cosmic rays. To be clear, define 
unit vectors in the directions of cosmic rays, Uj, AGNs, Vj, and Monte Carlo 
simulated cosmic rays, w^. A correlation si^na/ should then appear near 1.0 in 
the distribution of dot-products: Ui-Vj (if magnetic field deflections are mod- 
est). The index "i" runs over the cosmic rays in the data sample. For each 
value of "i", only the AGN catalog source (index "j") giving the maximum 
value of: u^Vj contributes to the distribution ^ . The simulated distribution 
of random background comes from the analogous distribution of: w^-Vj where 

^ Thus each cosmic ray has one entry in the c/oi-product distribution. This choice is 
consistent with each cosmic ray having one source. As only the AGN source nearest 
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index "/c" now runs over the sample of Monte Carlo simulated cosmic rays. As 
with the cosmic ray events, only the AGN catalog source (index "j") giving 
the maximum value of: w^-v, contributes to the distribution. 



3 Cosmic Ray and AGN Selection 



A few choices have been made in the comparison of AG AS A data and catalog 
of AGNs from RXTE. These are described here. 
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Fig. 1. The curves, from a GZK model [28], show the predicted fraction of cosmic ray 
events with source redshift, z < z-max versus cosmic ray energy for a selection of z-max 
distances: 0.003 < Zj^ax < 0.08. The GZK model assumed proton primaries with 
a power law spectrum at the source oc No enhancement factor of increasing 

source density with redshift, z, was included. The result was derived from Fig.6 of 
Ref. [28]. 



The AGASA data have energies, E > 40EeV and populate values of declina- 
tion: — 10° < Dec < 80°. As noted above, the steep cosmic rays spectrum, 
oc E~^, and modest number of events: 57 with E > 40EeV and 29 (just over 
half) with E > 53EeV led us to consider three (overlapping) bins in energy: 
E > iOEeV, E > 53EeV and E > lOOEeV. The last was to see if there are any 
correlations with the AGASA super-GZK events. Except for the E > lOOEeV 



in angle to the cosmic ray is chosen this can result in possible misidentification in 
the case of large source density. 
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selection, most of the cosmic rays are predicted, at least under the assump- 
tion of proton primaries [28], to originate at values of redshift, z > 0.01 ^ ; see 
Fig.l. 



RXTE efficiency estimated from Barger et al [astro-pfi/0410527] 
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Fig. 2. The curves show the estimated RXTE source detection efficiency, i.e. the 
fraction of RXTE AGNs [26] to all AGNs[30], VS source distance resulting from 
the RXTE instrument detection threshold. As noted in the text, two categories 
(or definitions) of AGNs are considered: o//-AGNs with RXTE 3-20keV intrinsic 
luminosities, -L3-20 > 10^^ ergs/s shown as "x", and broadline- AGNs with 2-8keV 
intrinsic luminosities, L2-8 > 10^^ ergs/s shown as "+". For redshift (distances) 
with full RXTE source detection efficiency then the fraction is 1. Parenthetically, 
we would obtain lower efficiencies (for the a//-AGN category) if we were to use the 
AGN number density VS X-ray luminosity deduced by the RXTE experiment [26] 
and shown as "*" . 

To match the AGASA acceptance, we selected AGNs with -10° < Dec < 80°. 
We have also made selections on the redshift of the AGNs to consider sources 
only with RXTE source detection efficiency ^50%. 

The estimate of the RXTE source detection efficiency involves two issues: 

(1) the RXTE instrument source detection threshold {i.e. the selection bias) 
VS redshift from Fig.l of Ref. [26], 

(2) the number density of AGNs VS redshift and intrinsic X-ray luminosity 
from Table 2 of Ref. [30]. 

Motivated by Ref. [30] , we divide the AGNs into the two categories: a//- AGNs 

^ This corresponds to a distance r 42 Mpc 
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and broadline-AG'Ns. For the all-AGN category wc require that the X-ray 3- 
20keV intrinsic luminosity ^, L3_2o ^ 10"^^ ergs/s, to match the RXTE data. 
With this intrinsic luminosity threshold the estimated a//-AGN number den- 
sity is 4.2 X 10~^ Mpc~^ consistent with the RXTE source density determina- 
tion of ~ 5 X 10~^ Mpc~^ [26] . For the broadline- AGN category we require that 
the X-ray 2-8keV intrinsic himinosity, L2-S > 10^^ ergs/s as this selects X-ray 
sources that are likely to be AGNs based purely on energetic grounds [29]. 
With this intrinsic luminosity threshold the estimated broadline-AGN number 
density is ~ 2 x 10~^ Mpc~^. 

Combining the RXTE detection threshold with our definition of two categories 
of AGN (above), we obtain the fraction of each AGN category VS redshift. 
This is shown in Fig. 2. Based on this result we restrict the redshifts for the all- 
AGN category to z < 0.005 and the redshifts for the broadline- AGN category 
to z< 0.03. 



4 Cosmic Ray— AGN Comparisons 

Plots of the distribution of rfoi- products (see definition in text) for the a^/-AGN 
selection are shown in Fig. 3. A plot of the AGASA cosmic ray and RXTE 
AGN directions are given in Fig. 4. The analogous plots for the broadline- 
AGN selection are shown in Fig. 5 and 6. A comparison of Figs. 4 and 6, 
shows two events shared between the two selections. Independently we have 
verified that all RXTE AGNs with redshift z < 0.03 satisfy at least one of our 
two AGN categories. 

The plots of the dot-products for the a//- AGN selection, Fig. 3, shows a small 
excess in the bin nearest to 1 for the AGASA event selections E > 40EeV, 
and E > 53EeV. For this bin (i.e. dot-product > 0.975) the excesses are ~ 1.1 
and ~ 1.7 standard deviations for the AGASA event selections E > 40EeV, 
and E > 53EeV respectively. If this correlation is valid, then it could provide 
experimental information to bound the magnetic deflections of extra-galactic 
cosmic rays. 

To see if the a//-AGN category excesses are consistent with e.g. typical GZK 
models, see Fig. 1, we estimate the RXTE AGN catalog efficiency as follows: 

(1) 90% sky coverage of the ~ 83% of the sky surveyed ^ ; 

^ For our study we relate the RXTE 3-20keV intrinsic luminosities, in ergs/s, 

to 2-8keV intrinsic himinositics, -L2-8 in ergs/s using: L2-8 ~ -^3-20 /2; private 
communication from Sergey Sazonov. 

^ This corresponds to the sky fraction outside a 10° avoidance zone about the 
galactic plane 
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Fig. 3. The plots show the distribution of doi-products (see definition in the text) 
for the aZ/-AGN selection: (top left) with cosmic ray energies E > 40EeV, (top 
right) with cosmic ray energies E > 53EeV, and (bottom) with cosmic ray energies 
E > lOOEeV. The curve on each figure shows the Monte Carlo random background 
normalized to the number of entries in each plot. 

(2) 87% estimated completeness factor; 

(3) ~ 77% estimated average all- AGN source detection efficiency (from Fig. 2). 

This yields an lower bound estimate for the a//-AGN RXTE efficiency of 
~ 50%. However for the all- AGN category (only 5 sources, see Fig. 4) it is 
likely that the global (redshift independent) RXTE efficiency factors are not 
appropriate ^ . To obtain an upper bound estimate for the all-AGN RXTE ef- 

^ In particular assuming an average AGN source density of 4.2 x 10^"^ Mpc^^ (see 
above) and the RXTE global efficiency [26] of 0.9 x 0.83 x 0.7 « 52%, the predicted 
number of nearby {z < 0.005) RXTE AGNs is approximated half those observed. 
While the small number of AGNs makes this weak statistically, it is nevertheless 
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Fig. 4. The figure shows the map of RA — Dec for the AGASA data and the 
AGNs from the a//-AGN selection. The AGASA data are plotted in hlue{*) for 
40EeV< E <53EeV, green{w) for 53EeV< E <100EeV, and red{U) for 100EeV< E. 
The RXTE AGNs are plotted as hlack{o). The galactic plane is drawn as a dotted 
line. 

ficiency we assume the global RXTE efficiency is ~ 100%. Then the estimated 
(upper bound) all-AG^ RXTE efficiency is ~ 77%. 

The estimated number of cosmic ray:a//-AGN coincidences is: the number of 
cosmic rays (with sources in a given redshift region) times the average a//-AGN 
source detection efficiency (for the same redshift region). Thus the estimated 
number of cosmic rays from the z < 0.005 region is obtained by dividing the 
excess counts, Fig. 3, by the a//-AGN efficiencies to obtain: 2.8/0.77 ~ 2.8/0.50 
or 3.6 ~ 5.6 events and 4.0/0.77 ~ 4.0/0.50 or 5.2 ~ 8.0 events respectively. As 
fractions of all the observed cosmic rays these are: 3.6 ~ 5.6/57 (or 6.3 ~ 9.8%) 
and 5.2 ~ 8.0/29 (or 18 ~ 28%) respectively. These fractions are somewhat, to 
significantly (depending on the a//-AGN RXTE efficiency), in excess of typical 
GZK models assuming proton primaries, see Fig. 1. 

Finally we note that the AGASA/HiRes cosmic ray quartet (or possibly quin- 

consistent with: a RXTE global efficiency of ~ 100% for nearby {z < 0.005) AGNs 
and/or with a local over-density of AGNs. Although these estimates were based 
on the AGN number density VS X-ray luminosity from Ref. [30] the AGN number 
density VS X-ray luminosity deduced by the RXTE experiment [26] gave a similar 
result. 
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Fig. 5. The plots show the distribution of doi-products for the hroadline-AGH 
selection: (top left) with cosmic ray energies E > 40EeV, (top right) with cosmic 
ray energies E > 53EeV, and (bottom) with cosmic ray energies E > lOOEeV. The 
curve on each figure shows the Monte Carlo random background normalized to the 
number of entries in each plot. 

tet[31]) cluster, at RA ^ 169.1°, Dec ^ 56.3° [15], is near one of the the 
RXTE AGNs at: RA 179.3°, Dec 55.23° and redshift z = 0.0035. In contrast, 
there is no close correlation in the all-AGN selection with any of the AGASA 
super-GZK events [4] plotted in Fig. 4 (more below). 

The plots of the (iot-products for the broadline-AGN selection. Fig. 5, are con- 
sistent with random background. If we assume the cosmic rays are primarily 
protons, then we can use a model such as Fig. 1 to estimate the number of 
cosmic rays expected from sources with redshifts z < 0.03. Then e.g. for the 
AGASA selection E > 53EeV, we expect ~ 60% to originate from sources 
with redshifts z < 0.03 or ~ 17.4 events. However the number that should 
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Fig. 6. The figure shows the map of RA — Dec for the AGASA data and the AGNs 
from the broadline- AGN selection. The AGASA data are plotted in blue{*) for 
40EeV< E <53EeV, green{w) for 53EeV< E <100EeV, and red{U) for 100EeV< E. 
The RXTE AGNs are plotted as hlack{o). The galactic plane is drawn as a dotted 
line. 

appear in dot-product bins near 1 depends on the RXTE AGN catalog effi- 
ciency. Similar to the evaluation above, we obtain an overall broadline-AGN 
RXTE efficiency of ~ 54%. Thus we should observe a signal as an excess of 
~ 9.4 events. Unfortunately in absence of a signal signature {i.e. dot-product 
bins in excess of random background) or of a bound on magnetic field deflec- 
tions, any statement on lack of excess depends on the assumed dot-product 
range. That said, assuming any signal would appear at dot-products > 0.95 
then ~ 9.4 events should result in a ~ 1.9 standard deviation excess. With 
the AGASA statistics and our current knowledge of cosmic ray deflections by 
magnetic fields, no strong conclusion can be drawn. 

The final issue is the evidence for, or against, correlations between the RXTE 
catalog of AGNs and the most energetic AGASA events. To investigate this, 
we show in Fig. 7 all of the AGNs from the RXTE catalog with z < 0.03 
and all of the AGASA events above lOOEeV as updated on the AGASA web 
site [32]. Now with 11 super-GZK events: 3 have dot-products > 0.975, ~ 3 are 
close to the galactic plane (region unobserved by RXTE) and the remaining 
5 do not correlate well {e.g. dot-product ^0.95) with the RXTE catalog of 
AGNs. 
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Fig. 7. The figure shows the map of RA — Dec for the updated hst [32] AGASA 
events with energies E > lOOEeV and the RXTE catalog of AGNs with redshift 
z < 0.03 without restriction on AGN dechnation. The AGASA events are plotted 
in red{M). The RXTE AGNs are plotted as black(o). The galactic plane is drawn as 
a dotted line. 

While the 3 AGASA super-GZK events that are close to RXTE catalog AGNs 
are consistent with random background, on inspection these events all have 
(iot-products > 0.99. In this case the expected random background is ~ 1.6. 
Furthermore one of these events, with energy E = 122EeV and RA 176.0°, 
Dec 36.3°, is close to the group of very nearby {z < 0.005 in Fig. 4) AGNs. 
The closest correlation is with the RXTE AGN at: RA 182.7°, Dec 39.45° and 
redshift z = 0.0033. 

Of the 5 AGASA super-GZK events that do not correlate well with the RXTE 
AGNs, 4 lie far from the galactic plane and have energies well above lOOEeV. 
Thus for proton primaries, based on Fig. 1 these should originate at redshifts 
^^0.01. For the broadline-AGN category the RXTE average source detection 
efficiency is then ~ 100% based on Fig. 2. Thus for these AGASA events we 
expect ~ 4 X 0.9 X 0.87 = 3.1 correlations with the RXTE catalog of AGNs 
whereas we observe zero. Furthermore, the Poisson probably of then observing 
zero is small: 4.4%. 

In contrast if the a//-AGN category is the more appropriate source of super- 
GZK events, then the RXTE average source detection efficiency based on Fig. 2 
is significantly less than 100%, particularly for source redshifts to z^O.Ol. In 
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this case, the a//-AGN category the RXTE average source detection efficiency 
is estimated at ~ 53%, resulting in an overall RXTE catalog efficiency of 
~ 34%. Thus we expect approximately 11 x 0.34 = 3.7 correlations (with 
sources to z^Q.Ql). The Poisson probability to observe < 1 correlation (one 
correlation was observed with the RXTE AGNs to 2<0.01) is 11.6% ^ However 
if the AGASA energies are overestimated (with respect to the energy scale of 
Fig. 1) then some of the AGASA events with energies closest to lOOEeV could 
originate at redshifts z > 0.01. If we extend the possible RXTE AGNs to red- 
shifts of 2;;^0.02 then the RXTE average source detection efficiency decreases 
to ~ 33% (because most of the lower X-ray luminosity AGNs are unobserved 
by RXTE) resulting in an overall RXTE catalog efficiency ~ 21%. Thus we 
expect 11 X 0.21 = 2.3 correlations and we observe two. The (new) additional 
correlation is between an AGASA event with E = 120EeV and a RXTE AGN 
at redshift z — 0.016. The Poisson probabihty to observe two correlations is 
~ 27%. 



5 Summciry 

We have searched for correlations between the published list of the highest 
energy events from the AGASA experiment [4,32] and the RXTE catalog of 
AGNs [26]. Two categories of RXTE AGNs were considered: all-AGNs with 
RXTE 3-20keV intrinsic luminosities, ^ 10^^ ergs/s, and broadline- 

AGNs with 2-8keV intrinsic luminosities, 1^2-8 ^ 10^^ ergs/s motivated by 
the analysis of AGN evolution in Ref. [30]. To retain RXTE source detection 
efficiencies ^50%, source redshifts of 2; < 0.005 and z < 0.03 were required for 
the a//-AGN and broadline- AGN categories respectively. 

No correlations were observed between the AGASA events and the broadline- 
AGN category of RXTE AGNs even though this category of AGN is most 
luminous in X-rays and even though the source density for this category of 
AGN is favored by some analyses [12,13] of the highest energy cosmic rays. 

In contrast, possible correlations were observed between AGASA events and 
the most inclusive, aZ/-AGN, category of RXTE AGNs. We note that while 
not statistically conclusive, one of the nearby RXTE AGNs correlates with 
the AGASA/HiRes quartet event cluster [15] and one correlates with one of 

^ If we also include the 320EeV event from the Fly's Eye [33] then we expect 12 x 
0.34 = 4.08 correlations (with RXTE AGNs to z^O.Ol) and we observe one. Now 
the Poisson probability to observe < 1 is 8.6%. Anecdotally the Fly's Eye event is 
very close, ~ 3.0°, to one of the RXTE sources at RA 88.8° Dec 46.3° and redshift 
z=0.02. If this is a true correlation, then the proton nature of the cosmic ray and/or 
the measured energy of the cosmic ray are in question. 
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the AGASA super-GZK events [32]. 



Additional data would help confirm, or refute, the interesting possibility of 
highest energy cosmic ray-AGN correlations. 
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